Effect of an inhomogeneous exchange field on the proximity effect in disordered 
superconductor-ferromagnet hybrid structures 
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We investigate the effect of an inhomogeneous exchange fieid on the proximity effect in 
superconductor-ferromagnet hybrid structures within the quasiciassical theory of superconductivity. 
As an example we study a superconductor-ferromagnet bilayer with an in-plane spiral magnetic 
order in the ferromagnet. The superconducting proximity effect induces in this case triplet corre- 
lations in the bilayer that are sensitive to the local quantization axis of the exchange field in the 
ferromagnet. The coexistence of singlet and triplet pair correlations in the bilayer results into a sen- 
sitivity of the superconducting transition temperature on the spatial variation of the exchange field 
in the ferromagnetic layer. We show that the inhomogeneity also tends to suppress the oscillating 
behavior of the pair amplitudes in the ferromagnet. 

PACS numbers: 74.45.+C, 74.62.-c, 74.78.-w 



I. INTRODUCTION 

The coexistence of superconductivity and magnetism 
is a long-standing issue and has gained recently a lot 
of attention due to experimental progress. One exam- 
ple is the discovery of ferromagnetic superconductors^. 
Typically, magnetism tends to suppress singlet super- 
conductivity via two mechanisms of pair-breaking : (i) 
the orbital pair-breaking effect by action on the electron 
charges; (ii) the paramagnetic pair-breaking effect by ac- 
tion on the electron spin via the Zeeman coupling. 

In contrast to the above-mentioned example of fer- 
romagnetic superconductors, in hybrid supcrconduc- 
tor/fcrromagnct (S/F) structures there is in general no 
coexistence of the magnetic and superconducting long 
range orders (assuming the ferromagnetic exchange field 
vanishes in S and the pairing interaction is repulsive or 
vanishingly small in F). Nevertheless, the influence of the 
magnetism on the superconductivity manifests itself near 
the S /F interface through the superconducting proximity 
effect: Cooper pairs can penetrate a certain distance into 
the ferromagnetic material. As a result of the exchange 
splitting of the Fermi surface in the ferromagnet, the 
Cooper pairs acquire a finite momentum and this causes 
spatial oscillations of the pair wave function in the F part. 
The oscillation of the pair amplitude in F shares many 
similarities^ in its origin with the Fuldc-Fcrrcll-Larkin- 
Ovchinnikov oscillations^ of the superconducting order 
parameter predicted to occur in the systems where the 
paramagnetic pair-breaking mechanism is dominant. In 
addition, the pair-breaking effect of the exchange field 
results in a much shorter penetration of singlet super- 
conducting pair correlations in the ferromagnet than in 
a nonmagnetic metal. 

In S/F bilayer systems a non-monotonic behavior is 
met again in the dependence of the superconducting 
critical temperature&i and in the dependence of the 
Josephson critical curren^SiSiifl on the ferromagnet thick- 
ness df or on the amplitude of the exchange field (for 
a recent review see Ref. ^] and references therein). 



Most of the quantitative investigations of these non- 
monotonic behaviors are based on the quasiciassical the- 
ory of superconductivity^^ which provides the simplest 
framework to study the inhomogeneity of pair correla- 
tions near the S/F interface. In the vaste majority of the 
theoretical works the exchange field J in the F layer is 
considered homogeneous for the sake of simplicity. How- 
ever in reality often the magnetic system is character- 
ized by the presence of an inhomogeneous magnetiza- 
tion, leading to a domain structure. In this paper we 
address specifically the influence of such domain walls on 
the proximity effect. 

Recently, Rusanov et alH have observed in S/F bi- 
layers a quantitative dependence of the superconducting 
critical temperature T c on the domain state of the ferro- 
magnet. In the presence of domain walls, T c is found to 
be enhanced compared to the T c obtained in the absence 
of domain walls. As the domain walls in Ref. Q are 
argued to be of the Neel type with an in plane magnetic 
moment, the orbital pair- breaking effect is here negligi- 
ble. The importance of the domain walls for the proxim- 
ity effect has also been pointed out in Ref. Il5i 

The above-mentioned effect is reminiscent to a similar 
effect in magnetic superconductors with an inhomoge- 
neous exchange field (for a review see Ref. E and ref- 
erences therein). As was first hypothesized by Matthias 
and SuhUi, one expects that the Cooper pairs experience 
an exchange field averaged over the superconducting co- 
herence length (£ s ) scale, which leads effectively to a re- 
duced pair-breaking effect for domain wall sizes compara- 
ble to or smaller than £ s . It is however worth mentioning 
that this qualitative picture is borrowed from the physics 
of coexistence of magnetism and superconductivity, and 
its applicability stricto sensu to hybrid structures may be 
questioned. One goal of this paper is to point out and 
investigate in depth the mechanism responsible for the 
sensitivity of the (singlet) superconductivity on the di- 
rectional changes of the exchange field in the S/F hybrid 
structures. 

The consideration of a full domain structure, i.e. with 
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FIG. 1: The exchange field rotates in the ferromagnetic layer 
in the plane (xy) with a constant wavevector Q. The period 
of rotation in the y direction is 2tv/Q. 



an alternation of regions with fixed magnetizations and 
domain walls represents a formidable theoretical task. 
This calls naturally for the study of simplified models 
in a preliminary step. One possibility is to study a local 
inhomogeneity of J in the vicinity of the S/F interface. It 
has been foundi2ii2i2S that superconducting triplet corre- 
lations with an unusually long penetration length in the 
ferromagnet arise through the S/F proximity effect. Such 
long range triplet correlations have been shown also to 
be characteristic for F/S/F trilayer structures with non- 
collincar F momentsi 21 ! 22 It was concluded 18 ' 19 i 20 i 21 i 22 
from these two models that the sufficient ingredient for 
the existence of the long-range triplet components is a 
change in direction of the F moment. 

In a recent paper j2£ we have studied an S/F bilayer 
within a model of an in-plane rotating magnetization 
(spiral order—) in the F layer (see Fig. QJ. This model 
is consistent with the fact that the domain structure 
is expected a priori to appear in the layer rather than 
across the layer. We founci 2 ^ that the long-range triplet 
components are not induced, although the exchange field 
is inhomogeneous. To understand this peculiarity, it is 
worth noting that the models of moment inhomogencities 
studied in the Ref. led to one-dimensional 

spatial dependences of the pair correlations: the inho- 
mogeneity of the moment was always considered across 
the layers, i.e. competing with the inhomogeneity due to 
the proximity effect itself. Contrary to this, the moment 
inhomogeneity within the in-plane spiral order model oc- 
curs in the transverse direction, so that the problem is 
intrinsically two-dimensional. 

As pointed out in Ref. [2f|, the anomalous Green func- 
tion / and the Green function g, which are 2x2 
matrices in spin space, acquire both a nontrivial spin 
structure through the proximity effect. It has been 
realized22iSi2Si2i that this always results in short-range 
triplet components (i.e. even in absence of inhomogene- 
ity of J) induced together with the singlet component in 
the ferromagnet near the S/F interface. For example, the 
difference in the T c in F/S/F trilayers in the parallel and 
antiparallcl configurations of the F moments 2 ^ stems pre- 



cisely from these short-range triplet components— The 
same mechanism is responsible for the sensitivity of T c on 
the degree of inhomogeneity in the model of the rotating 
magnetization— 

In the next section of this paper, we develop our gen- 
eral approach within the framework of the quasiclassical 
theory of superconductivity. We derive the conditions re- 
quired for the presence of long-range triplet components. 
In Sec. Ill we investigate how the S/F proximity effect 
depends quantitatively on the degree of moment inhomo- 
geneity in an S /F bilayer structure within the spiral order 
model. We provide the details of the calculations leading 
to the results already presented in our short paper—. In 
addition to Ref. |23, we present the results concerning 
the spatial dependences of the singlet and short-range 
triplet components and discuss quantitatively the condi- 
tions which may increase or reduce the triplet amplitude. 



II. SPIN STRUCTURE OF THE GREEN 
FUNCTIONS INDUCED BY THE S/F 
PROXIMITY EFFECT 

In this section, we show that it is possible to capture 
general features of the S/F proximity effect independent 
of the specific geometry of the system under consider- 
ation. For example, a spin splitting in the local den- 
sity of states 2 ^ is straightforwardly found to result from 
the production of the triplet components near the S/F 
interface— This spin imbalance is obviously accompa- 
nied by the penetration of a spin magnetization in the 
superconductor—. We shall also put forward the condi- 
tions and the physical reasons for the production of the 
long-range triplet components. 



A. Equations in Nambu-Gor'kov space 

We discuss the S/F proximity effect within the frame- 
work of the quasiclassical theory of superconductivity 
within the Nambu-Gor'kov formalism. The particle- 
hole quantum-mechanical coherence is expressed by a 
2x2 matrix structure of the Green function in Nambu- 
Gor'kov space, 
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(1) 



The functions g and / are the particle-hole conjugates 
of the Green functions g and /. We assume that the 
S/F bilayer we consider is in the diffusive limit, in which 
case the Green function g(R,ui n ) depends on a spatial 
coordinate R, and on the Matsubara frequencies u> n = 
7rT(2n+I ) (n an integer, T the temperature). The Green 
function obeys the Usadel transport equation for diffusive 
systems^ 



i lu„t 3 - A - J(R) • <7, g 



D; 



■Vi (ffV 3 -5) = (2) 
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with the normalization condition 
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(3) 



where is the diffusion constant tensor, fn Eq. J2J), 
we used the Einstein convention for summation over re- 
peated indices = x, y, z). We assume in our cal- 
culations for simplicity isotropic diffusion tensors, = 
DSij. The diffusion constant D in the superconductor 
differs a priori from the diffusion constant in the fer- 
romagnet (an index characterizing both values is intro- 
duced below). We assume the same diffusion constants 
for both spin projections in the ferromagnet. The su- 
perconducting order parameter A and the spin matrix a 
have the form 



A 



A 
A 



and <t 



a 

a* 



(4) 



Here <x is the vector of spin Pauli matrices. We consider 
a superconductor with conventional s-wave pairing, so 
that A = A*. In the ferromagnet, the superconducting 
order parameter vanishes, A = 0. The proximity effect 
manifests itself in a nonzero pair amplitude / ^ in 
the ferromagnet. In the superconductor, the exchange 
field vanishes, J (R) = 0, expressing the non-coexistence 
of superconducting and ferromagnetic orders. A non- 
zero magnetization <5M(R) in the superconductor may 
be induced by the inverse proximity effect. 

The Nambu-Gor'kov representation includes some re- 
dundancy which manifests itself in the fundamental sym- 
metry relations^ g(R,u) n ) = <?*(R, cj„), /(R, u) n ) = 
f*(R,u> n ), and 



g(R,-u n ) = g^(R,u n ) 
f(R,-u n ) = -/ tr (R,cj„). 



(5) 
(6) 



Here and g tr are respectively the adjoint and trans- 
posed matrices. 



B. Spin structure 

The Green function g and the anomalous Green func- 
tion / are 2x2 spin matrices. For example, the anoma- 
lous Green function is written as 



f _ ( /tt fn 
\ /it fu 



(I) 



The order parameter in spin space for singlet pairing 
reads 



A = A s ia v 



(8) 



Generally 


30 


31 


32 


33) 



/ = [fa + (ff ' o)\ io y = 



-ftx + ifty ftz + fa 
ftz fs ftx ~t~ ifty 



(9) 



where f t = (ftx, fty, ftz) is the triplet pairing vector and 
f s is the singlet pair amplitude. J defines the spin- 
quantization axis (while in Sec. II we choose to coincide 
with the z-axis). Thus, the probability to find a pair in 
the triplet state with zero spin projection on the quanti- 
sation axis z is proportional to \ftz\ 2 ■ When the vector f t 
is directed along the exchange field J (i.e. ft y = ftx — 0), 
the Cooper pair spin is perpendicular to J and does not 
contribute to the spin paramagnetic susceptibility. When 
the triplet vector t t is noncollincar with J, it means that 
the amplitudes ±f tx + ift y of the states with spin projec- 
tions ±1 on the quantization axis are non-zero. In this 
case, the Cooper pairs contribute to the susceptibility. 
As for the Green function g, we adopt the notation 



1 g x + %g v go - g z 



(10) 



with g = (g x , g y , g z ) the vector part of g, and go its scalar 
part. Also, we expand the particle-hole conjugates in the 
spin space as 

f = i<r v f a -(ff<r) , and g = g - <r y g ■ aa y . (11) 



C. Normalization condition in spin-space 

The normalization condition |J5J written for the differ- 
ent components of the Nambu space yields the relations 



gf + fg = o 
g 2 + ff = -ttVq 



(12) 
(13) 



where oo is the unit matrix in the spin-space. In equilib- 
rium, it can be shown from Eqs. I|12|) and (|13fl that for 
physical solutions necessarily the condition 



Tr (<?) = 0, i-e. go 



-9o 



(14) 



holds, which expresses the particle-hole symmetry in the 
hybrid S/F structure. 

Using the fact that the matrices (<7o, <t x , cr y , a z ) form a 
basis for the spin matrices, we find from Eqs. I|12|) and 
(|14f> that the different components of / and g in spin 
space obey the following conditions 



ff(g-g) = 0, 
fs (g - g) = if* x (g + g) 

Eq. (C3J) yields 

gl+E 2 ~fafa+ ft -ft - 

25og + fsit - fsft + i ft x f t = 







(15) 
(16) 

(17) 
(18) 



with g 2 = g\ + g^ + g\ = g • g |g| if the components 
of g are complex numbers) . By combining the last equa- 
tion l|18|) with its particle-hole conjugate, we obtain the 
relations 



go(s -g) = i ft x ft 
5o(g + g) = / s ft - fsit- 



(19) 
(20) 
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As Eqs. (CSJ) and (TBJ follow from and (J20J), the 

normalization condition © together with l|14|l leads to 3 
independent equations (|17|) . 1)19(1 and H20|l in spin-space. 
Note that from Eq. (JTHJl and (J2UJ) follows g 2 = g 2 . How- 
ever, in general g ^ g. According to Eq. (|19|l . the equal- 
ity holds when f t || f*. 

Moreover, we see that necessarily g = g = when 
t t = 0. This implies a spin-independent density of states 
in the ferromagnet in the absence of the proximity ef- 
fect. This finding reflects the fact, that changes in the 
normal-state density of states as a result of an exchange 
splitting are small in the expansion parameters of quasi- 
classical theory as long as the exchange splitting is small 
compared to the conduction band widths for both spin 
directions. Concordantly, we assume spin-independent 
diffusion constants for consistency. In the opposite case, 
when the splitting is large compared to all low-energy 
scales of the problem (T, A s ), the spin splitting must be 
taken into account as spin-dependent dispersions before 
applying the quasiclassical approximation (see Ref. l34l) . 

D. Quasi-classical transport equations in spin-space 

We now derive the transport equations for the spin 
components of / and g. For the component / the equa- 
tion J2J) yields 

— Vj (gVjf + fV 3 ~g) + 2i to n f = Ag - gA + 

+J-0-/-/J-0-*. (21) 

The system of equations in the spin space reads (taking 
into account that go — —go) 

-Vj [goVjfs - fsVjgo + g • Vjf t - f t ■ Vjg] 

+2i uu n f s = -2g A s + 2 J • f t (22) 

— Vj [ffoVjf t - ftVjgo + gVj.f s - /sVjg + igx Vjft 

TT 

+i Vjg x f t ] + 2i u> n i t = -A s (g + g) + 2 J/ 8 . (23) 

It is supplemented by the self-consistent equation 
defining the s-wavc order parameter A from f a which 
reads in the dirty limit for a weak-coupling 

A S (R) = AT^/ S (R, W „) (24) 

n 

with A the pairing interaction constant. Using the equa- 
tion 

^En = l + 4 

„ Kl A T 

relating the critical temperature T c o of the superconduc- 
tor without the proximity of a ferromagnetic layer and 
A, we can rewrite the Eq. I)24|l as 




Near the critical temperature T c , the pair amplitudes 
f 8 and f t are small and the Green function g deviates 
only slightly from its value (g = — wrco sgn LJ n ) in the 
normal state, so that the Usadel equations (|22|l - 123|l can 
be linearized and take the simpler formS 

(DV 2 -2\uj n \) f s = -27rA s + 2isgnK)J-f t (26) 
(DV 2 -2\u; n \) f t = 2isgnK)J/ s . (27) 



E. General features of the S/F proximity effect 

Within a particular geometry for the S/F hybrid struc- 
tures, we have to complement the transport equations 
with boundary conditions. Nevertheless, many general 
features of the S/F proximity effect such as the number 
and the nature of the nonzero spin components for / and 
g can be determined independently from the geometric 
effects of the boundaries. 

From the transport Eq. IL'-ili , and even clearer from Eq. 
(|27|l . it follows that the triplet vector f t is necessarily 
non-zero if the singlet component f s penetrates in the 
ferromagnet. Moreover, the triplet vector tends to align 
with the exchange field J, indicating that the dominant 
triplet component has zero spin-projection. Thus, the 
singlet component and the triplet component with zero 
spin-projection coexist always in the ferromagnet near 
the S/F interface. This is physically expected since these 
two pair correlations are energetically equivalent in the 
ferromagnet with respect to their interaction with the 
exchange field. Both components are characterized by 
short-range penetration lengths in the ferromagnet. 

On the contrary, if the triplet vector f t is non-collincar 
with J, it means that triplet components with nonzero 
spin-projection on J are produced. Since these corre- 
spond to equal-spin pairing, they are not limited locally 
by the paramagnetic interaction with the local exchange 
field and may have long-range scales in the ferromag- 
net. A misalignment between the triplet vector f t and 
the moment J occurs in presence of sudden changes in 
orientation of J. The reason is that ft obeys a differen- 
tial equation and its variations in orientation have thus 
to be relatively smooth. 

From Eq. (|20() it is clear that the counterpart for 
the production of triplet components (f t ^ 0) for / is 
the presence of g ^ for g. As a direct consequence, 
the density of states for the up and down spin projec- 
tions differs^ Furthermore, in the presence of long-range 
triplet components, the Green function g contains also 
off-diagonal (spin-flip) terms. As a result of the spin 
splitting in the density of states generated by the S/F 
proximity effect, a spin magnetization <5M is induced near 
the S/F interface. This magnetization leakage has been 
investigated recently in the Ref. |53 within a model con- 
sidering a fixed exchange field. The spin magnetization 
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induced by the proximity effect is given by 3 ^ 3 ^ 

SM(R) = 2iVoT^g(R,uv). (28) 

n 

Since the triplet vector f t is also induced in the supercon- 
ductor near the S/F interface via an inverse proximity 
effect, the vector g characterizing the magnetic correla- 
tions penetrates also in the superconductor according to 
the relation (J2DJ. 

As we show in the following, the sum over Matsubara 
frequencies in the expression (|28|) is in general nonzero. 
Indeed, as noticed in the papersi&SiiS^, in the diffusive 
limit the triplet components have the property to be odd 
functions of the Matsubara frequencies u> n , while the sin- 
glet amplitude f s is an even function of uj n . These prop- 
erties follow directly from the Pauli principle, which leads 
to the relation ©. The different components of g have 
also symmetry properties with respect to u n . Combining 
()14|l with Q, we obtain that go is an odd function of ui n , 
and 

g(-w„) = gK). (29) 

We note that the relations (|19|I - 12U|) between the vec- 
tors g and f t are simplified when f t x f t = 0. This condi- 
tion corresponds to unitary triplet superconductivity 3 !. 
It can be seen from the transport equations H22fl - H23|) 
that if the gap A s can be chosen real (and thus the sin- 
glet amplitude f s is real), then the triplet vector ft is 
purely imaginary, i.e. ft = f t * = —ft (taking into account 
that go = —go, i.e. go is purely imaginary). As a result, 
one obtains the simpler relations 

g = g, Sog = fA- (30) 

These simplifications are found also if ft || J (i.e. when no 
long-range triplet components are induced). Combining 
the relations (|30|) and l|29|) . we conclude that the spin- 
vector part g is an even function of uj n under certain 
circumstances, which demonstrates that the induced spin 
magnetization <5M is non-zero. 

Therefore, although there is no coexistence of ferro- 
magnetic and superconducting orders, magnetic and su- 
perconducting correlations do coexist in the vicinity of 
both sides of the S/F interface. It is worth noting that 
near T c , the singlet amplitude f s and the triplet vector 
amplitude are small so that g and thus <5M appear to 
be second order terms (see Eq. fflfy ). Accordingly, one 
expects that the induced spin magnetization <5M pene- 
trating the superconductor, which is negligibly small near 
T c , increases significantly by reaching temperatures well 
below T c . 

III. S/F BILAYER WITH A ROTATING 
EXCHANGE FIELD 

In this section, we study the proximity effect in a S/F 
bilayer within the model of a rotating exchange field in 



F (see Fig. nj. We derive analytical expressions for the 
spatial dependences of the singlet and triplet amplitudes 
near the superconducting critical T c as a function of the 
spiral wavevector Q. We calculate also numerically the 
dependence of T c on Q. A part of the results has been 
already reported in our short paper—. We shall see here 
that the magnetization induced in the superconductor by 
the triplet components reflects the inhomogencity of the 
exchange field J in the ferromagnet. 

A. Model and boundary conditions 

In the F layer (defined by the plane (x,y)), J rotates 
with an angle varying along the direction y, i.e. 

J(y) = J (cos Qy, sin Qy, 0). (31) 

It is straightforward to see that in the present geometry 
the component f tz of the triplet vector is zero (because 
J never points in this direction). It is worth mentioning 
that here z does not correspond to the spin quantization 
axis. To identify physically the different components of 
the triplet vector, one has to express ft in the local basis 
(X, Y, Z) where Z is the direction of the local exchange 
field J. According to the symmetry relations (see the 
previous section), it is sufficient to consider the positive 
Matsubara frequencies. We want to determine the de- 
pendence of the critical temperature T c on the rotation 
wavevector Q, and for this purpose we have to solve the 
linearized Usadcl equations. The y dependence of the 
moment J is eliminated in the right-hand side of the Eq. 
(|26[) - <|27[) by considering the new components 

/+ = (-ft x +ifty)e lQ y (32) 
/- = (f tx +ifty)e- iQy . (33) 

The new system of equations to solve takes the form 

(DV 2 - 2uj n ) f s = -2ttA s + i J (/_ - /+) (34) 
(DV 2 =F 2iDQd y - DQ 2 - 2u; n ) f ± = T 2i Jf s . (35) 

The diffusion constants are D = D s in the S layer and 
D = D f in F. 

The components of the triplet vector in the plane (cc, y) 
are obtained from / + and /_ with the relations 

ft* = \ {f-e lQy - f+e- lQv ) , (36) 

fty = \ % {f-e iQ y + f+e-^y) . (37) 

In the present case, the gap amplitude A s and the singlet 
amplitude f s can be chosen real. Then /_ = /?, i.e. 

ftx = i (-Im /+ cos Qy + Re /+ sin Qy) (38) 
fty = -i (Re /+ cos Qy + Im /+ sin Qy) . (39) 

It is clear that the amplitudes of the triplet components 
are purely imaginary. 
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The system of equations <|34|) - (I35|) ^ s supplemented 
by boundary conditions. The boundary conditions at 
the S/F interface (provided that J <C £f with ep the 
Fermi energy) for the diffusive regime have been for- 
mulated by Kuprianov and Lukichev for a small barrier 
transparency^. The general boundary conditions have 
been derived by Nazaro\*2&. Near T c , they are formally 
equivalent and reduce to 

6 d z f) sc = 7 £/ 8 z f) F , 7 = P.Z./P& (40) 

where p s and pf are respectively the normal-state resis- 
tivities of the S and F metals (this boundary condition 
follows from the continuity of the current at the inter- 
face), and 

tnt a„f) F = f) sc - f) F , 76 = RbA/pftf (4i) 

with Rb the resistance of the S/F boundary, and A its 
area. Here £ s = y D s /2ttT c q is the superconducting co- 
herence length in S and £/ = \/Df /2irT c o is the super- 
conducting coherence length in F. At the outer surfaces 
of the F or S layers [z = —df and z = d s ), the current 
through the boundary has to vanish, i.e. 



d z f = 0. 



(42) 



It is important to note that the present boundary con- 
ditions H4U|) - (|42fl do not couple the different spin compo- 
nents of /. 



B. In the ferromagnet 

In the ferromagnetic layer, the singlet amplitude and 
the triplet components f± are coupled through 



(v 2 -Q„£7 2 )/ 8 = iC 2 (/_-/+) 
V 2 =F 2iQd v -Q 2 - 0„C7 2 ) /± = =F 2i C, 2 f s . 



(43) 
(44) 



where ti = y/Df/J and 0„ = (2n + 1)T/T c0 . 

Since the geometry is periodic in the y direction, 
the components of the superconducting condensate wave 
function / can be expanded into a Fourier series. Using 
the boundary condition at the outer surface (z = —df), 
the components of / are sought under the form 



fl(y,z) = 



E f {p) 

P— — OC 



fi cosh [kf p (z + df)} e 



ipQy 



(45) 



where I = s, ±. Substituting these expressions in the set 
of equations (|43|) - (|44|l leads to the following linear system 



*7 

2 



2 kl-Q 2 



k 2 Qp—i 



( 

Ap) 





(46) 



with k 2 = k 2 p — fl n ^f 2 and Q p = pQ. The eigenvalues 
k 2 p are determined from the condition of zero determi- 
nant for the 3x3 matrix. For p ^ the three eigenval- 
ues k 2 p - (j = 1,2 and 3) and the associated eigenvectors 
(fa,p,ji f-,p,ji f+.p.j) take a complicated analytical form 
that we shall not give here. We note that for p ^ 0, the 
three components of each eigenvectors fi lP j are non zero. 

For p = 0, the three eigenvalues can be easily derived 
and have a simple form 



K fOe 



0„C/ 2 + £2i£j 2 r)-e and k 2 f03 = Q n ^ 2 + Q 2 

(47) 



with r) e = y/l — rj 2 + ieij for 77 < 1, e = ±1, rj = (jQ 2 
For r\ > 1, we have rj £ = —i\/if — 1 + ier) . The two first 
eigenvalues (index e = + or cquivalently j = 1, and e = 
— or j = 2) which are related by a complex conjugation 
for r) < 1 correspond to a short penetration length (at 
least at small Q) of the order of £j -C £/ f° r strong 
ferromagnets (J ^> 27rT c0 ). On the contrary, the third 
eigenvalue (fc 2 03 ), which is independent of the exchange 
field amplitude J, is real and gives a much longer decay 
length of the order of for the pair amplitude in the 
ferromagnet. The corresponding eigenvectors have the 
form 




Ve 



and 




(48) 



In the limit Q — > 0, then rj — > and rj e — > 1, and we 
find again from (|47J) the known eigenvalues for a fixed 
exchange fieldS^. 

The general solution of the system (|43|l - 144() satisfying 
the outer boundary condition can be written as 

+00 3 

/»(». *) = E E cosh [ fc / P ^ z + eipQ " 

p—— 00 j—i 

(49) 

where the three coefficients a } have to be determined 
with the help of the boundary conditions at the S/F in- 
terface. 



C. In the superconductor 

1. Triplet components 

In the S layer, there is no coupling by the equations be- 
tween the singlet and triplet components. The solutions 
for f± satisfying the boundary condition at the outer sur- 
face (at z = d s ) and being periodical in the y direction 
are straightforwardly derived and have the form 

f± (y, z) = f± cosh [ k U z ~ ^)] eipQy ( 5 °) 



where 



Mnlis 2 + Q 



2 



(51) 
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and the coefficients f±' have to be determined with the 
boundary conditions at the S/F interface. 



2. Singlet component 

In the S layer, the equation for the singlet pair ampli- 
tude f s 



(D S V 2 - 2u n ) f s 



-2wA, 



(52) 



is coupled to the self-consistency equation (|25fl . Due to 
the periodic geometry in the y direction, f s and A s can 
also be expanded into Fourier series 



f s (y,z)= J2 /i P) W e 



ipQy 



p— — oo 



A s (y,z) = A i P) (z)e 



ipQy 



(53) 



(54) 



Then, the Fourier amplitude fs (z) obeys the differential 
equation 

{d 2 z -Ql- fi n C 2 ) = -2nA^/D s (55) 
while the gap amplitude A^ (z) is given by 



A(^)ln^ = 2.TW^-^M 

n>0 V 



(56) 



In general, the coupled equations ((55() - l(56() can n °t be 
solved analytically due to the self-consistence. 

It is clear at this point that there is no mixing between 
the different Fourier components of the superconduct- 
ing condensate function / neither in the ferromagnet nor 
in the superconductor near the critical temperature. It 
means that each single Fourier component p (which de- 
termines a particular y dependence) is a solution of the 
system of equations which satisfies the boundary con- 
ditions. For each p, we obtain a different gap equation 
H56p. i.e. a different critical temperature T c {p). The solu- 
tion realized physically is the one which gives the highest 
critical temperature (i.e. which is energetically most fa- 
vorable). 

For p ^ 0, the singlet amplitude f s and the order pa- 
rameter A s are inhomogencous along the direction y even 
far from the S/F interface. The solution with p = cor- 
responds to a singlet amplitude f s which depends only 
on the spatial variable z both in the F and S layers. In 
this case, the components f± are also independent of the 
coordinate y (characterizing the spatial inhomogeneity of 
J). The y dependence of the triplet vector is then sim- 
ply revealed by the relations (|3*5|l - ll3T|) . The influence of 
the inhomogeneity of the exchange field on the singlet 
amplitude f s in the superconductor occurs near T c only 
through the boundary conditions at the S/F interface. 



D. Pair amplitudes 



(p) 



The remaining step is to determine the coefficients al- 
and with the boundary conditions at the S/F in 
tcrface. Let us define the short-hand notation 5 
fi p \z = 0). The condition gTJ yields 



(p) 



3=1 
3 



(p) 



(p) 
i 



where I = ± and 

•Aj = cosh(k fpj d f ) + Jbkf P j€f smh{k fpj d f ) 



(57) 
(58) 

(59) 



Then, the condition 14U|) considered for the triplet com- 
ponents only, combined with the Eq. I|58(l gives 



-VP) c ~\iv) _ n 



(60) 



with 



A® = Af+Cf\kU s t^h{k\ p d s )) 



i(p) 



(p) (U 



(61) 



c j = lk fpj £f 8inh(kf pj d f ) 



(62) 



The equations 1(57(1 and l(60() lead to a 3 x 3 linear system 
for the coefficients (if - As a result, we obtain straight- 
forwardly for p ^ 



,(p) _ r(p) N 3 
l j - 6 ° F-N 



(63) 



with N = h+ x h , and where h; (I = ±) and F are 
vectors whose components labeled by j = 1, 2 and 3 are 
given by 



F — f A 



(p) 



(P) 



h — f A VP - 

lH,j — JLp,j^j,l 



(64) 



For p = 0, the absence of coupling between the singlet 
amplitude f s ,3 and the two triplet amplitudes /±.3 for the 
third eigenvector in (|48|l yields that the long-range con- 
tribution is absent (a^ = 0). By the same occasion, the 
quantity A^} is independent of the index I = ± (which 
is dropped out in the forthcoming expressions). Accord- 
ing to the form of the eigenvectors associated with the 
eigenvalues k^ e in Eq. (148(1 (from now on we drop the 
index indicating that p = 0), the two triplet components 
/+ and /_ in the ferromagnet are related as / + = — /_. 
Since in addition the symmetry property /+ = holds, 
we have necessarily Re / + = 0. Then using the Eq. 138(1 - 
(|39(l . we find that the corresponding triplet vector has 
the spatial dependence 



ftx(y, z) = i Im /_ (z) cos Qy 
fty(y,z) =«Im/-(z) s'mQy, 



(65) 
(66) 
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i.e. ff follows everywhere the direction of the inhomogc- 
neous exchange field J(y). Thus, only the triplet vector 
component f t z with zero spin-projection on the local mo- 
ment J (direction Z) exists. One obtains that 



ftz(z) = f-(z) 



e=± 



a e e cosh[k f e (z + df)] (67) 



is characterized by a short penetration length in the fcr- 
romagnet, as expected from the physical arguments pre- 
sented in the Section II. In the F layer, ftz exhibits in 
addition oscillations in space as long as the eigenvalues 
k 2 £ are complex. The two remaining nonzero amplitudes 
a+ and a_ in the ferromagnet (for convenience we label 
the amplitudes by e = ± rather than j = 1, 2) obey the 
2x2 system determined by the Eq. ljST|) and JBOJl, and 
are given by 



a + 



A 2 



So 



V+AiA 2 + rf-AiAi 
V+AxA 2 + r\-A% Ai ' 



(68) 
(69) 



For r\ < 1, we have a + = a*_, while for ry > 1 the co- 
efficcnts a± are purely imaginary. We find from the Eq. 
()58(l that the triplet component with zero spin projection 
has the following spatial dependence in the S layer 

„ , . , , cosh \kt(z — cL)l 

f tz {z) = {a+A, - a-A 2 ) A, , , J 70 

cosh(fctd s ) 

with the real wavevector k t = y/ £l n £,7 2 + Q 2 ■ Finally, it 
is straightforward to see that for p = 0, the spin-vector 
part g (an even function of lo u here) of the Green function 
g determined from the singlet and triplet components 
according to the relation i|30|) can be written near T c as 



-fs(z)ftz(z)3(y). 



(71) 



For p = 0, the function Wo can be written as 



Wn 



f] + C 1 A 2 +ri_C 2 A 1 
V+AiA 2 +r)-A 2 Ai 



(74) 



In the limit Q — ► 0, we recover from this expression the 
formula (12) of the Ref. Q obtained for a S/F bilayer 
with a constant exchange field J in the ferromagnetic 
layer. All the informations characterizing the proximity 
effect between the S and F layers are contained in this 
real function W p . 

Following the Ref. 0, the solution of the Eq. (|55|l sat- 
isfying the boundary conditions at the outer surface (Eq. 
|g2J for z — d s ) and at the S/F interface (Eq. l(72^) is 
expressed using the Green function of the equation. As 
a result, one can write the singlet pair amplitude as 

/j%)=7r [ d G(y,z)A^(y)dy. (75) 
Jo 

with the Green function G given by 



G(y,z) 



(fc^TTTeo)- 1 



sinh(fc s d s ) + (W p /k s £ s )cosh.(k s d s ) 

vi(z)v 2 {y), z <y 



< v 2 {z)vi{y), y <z 



where 



Vx{z) = cosh(k s z) + (W p /k4 s )smh(k s z), (77) 
v 2 (z) = cosh(fc s [z — d s ]), (78) 



and 



k s — \ QnC 



(79) 



Combining the equations (|56|l and (|T5|> - (|79|) one obtains 
a single equation to be solved numerically. 



From this equation, we note that the spin magnetization 
induced in the superconducting layer exhibits the same 
inhomogeneity as the exchange field J(y) in the trans- 
verse direction y. 



E. Superconducting critical temperature T c 

Using the last boundary condition l|4()|) considered for 
the singlet amplitude and the expression l|63|l , we derive a 
relation between the derivative of the singlet component 
fs (z) and its value 5^ at the interface as 



ts& Y (z = 0) = W p 6% 



(/<) 



with 



C-N 
F • N ' 



(72) 



(73) 



F. Numerical results and Discussion 

The dependence of T c on the spiral wavevector Q for 
the different harmonic solutions is plotted in Fig. (for 
definiteness, we took here similar parameters as in Fig. 2 
of Ref. H that is T c0 = 7 K, 7 = 0.15, J = 130 K, 
£ s = 8.9 nm and £/ = 7.6 ran). As shown in this figure, 
the harmonic p = yields the highest T c . We found the 
same result for a large range of input parameters (e.g. 
for bigger thicknesses d s and df, not shown here). The 
solutions with p / correspond to a superconducting 
state in the S layer which remains inhomogeneous far 
from the S/F interface (corresponding to a Fulde-Ferrell- 
Larkin-Ovchinnikov stated) . They are energetically less 
favorable than the solution with p = which gives a 
singlet pair amplitude that is inhomogeneous only near 
the interface due to the proximity effect (the fact that the 
singlet correlations arc independent of the coordinate y 



d f =2.5nm 
d =7nm 




FIG. 2: T c versus the spiral wavevector Q for different har- 
monic solutions. The highest critical temperature at Q 7^ 
is obtained for p = 0. Here 75 = 0.3, df = 2.5 nm and d s — 7 
nm. 



for a y-indcpcndent order parameter simply reflects their 
isotropy in the spin space: singlet pairs do not see the 
directional changes of the exchange field). 

Henceforth, we discuss only the physically relevant so- 
lution p = 0. The superconducting critical temperature 
T c clearly increases with Q as shown in the Fig.[21(scc also 
the Fig. 2 and Fig. 3 of Ref. I23T) . This enhancement is 
observed for a large range of parameters we have tested, 
and thus is likely to be general. The dependence of T c 
as a function of df is shown in Fig. [31 for two different 
values of the dimensionless parameter Q£ s (here we took 
J = 20T c0 , 6 = £/, d s = 2£ s , 7 = 0.15 and 7b = 0). An 
important characteristic feature revealed by this figure is 
the evolution of the nonmonotonic dependence of T c (df) 
in favor of a monotonic one in the presence of a moment 
inhomogeneity in the ferromagnet. 

As an exacerbation of this nonmonotonic behavior of 
T c (df), it is well-knowni that within a particular choice 
of parameters the superconducting critical temperature 
may even jump to zero for a finite range of thicknesses 
df with a reentrance of superconductivity for higher val- 
ues of df. We investigated the influence of the moment 
inhomogeneity within this particular r ang e of df where 
T c = 0. As shown in the Fig. 2 of Ref. |22, we found the 
restoration of the superconductivity above a threshold 
value for the spiral wavevector Q. It is worth mentioning 
that a similar reentrance of superconductivity caused by 
a moment inhomogeneity is known in some rare-earth in- 
tcrmetallic compounds where a magnetic order competes 
with a superconducting order (see Ref. Il6|) . In these com- 
pounds superconductivity is accompanied by a transition 
from the homogeneous ferromagnetic state to the crypto- 
ferromagnetic stated. However, it is important to note 
that this similar reentrance behavior in bulk systems and 
in S/F hybrid structures does not result from a similar 
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FIG. 3: T c versus d f for Q£ s = and Q£ s = 1. The 
characteristic non- monotonic dependence of T c (df) found at 
Q = is suppressed by the moment inhomogeneity in F. Here, 
J = 20T c0 , £ s = £/> d s = 2£ 3 , 7 = 0.15 and j b = 0. 



physical mechanism. In S/F structures this behavior re- 
sults from a non-local influence of the magnetic order 
on the superconducting order (since both orders are as- 
sumed to be spatially separate) based on the Andrcev 
reflection process. According to the Sect. II, this process 
is characterized by the coexistence of singlet and triplet 
pair correlations near the S/F interface. 

In the Fig. 0] we have represented the dependences 
T C (Q) for different values of the parameter 7^ (the other 
parameters are the same as in Fig. 2 and Fig. 3 of 
Ref. [2^). The threshold value for Q corresponding to 
the switch from the normal state to the superconducting 
state depends as expected on the quality of the S/F in- 




Yb =0.4 



/Yb=0-2 



Yb =0.0 



12 3 

FIG. 4: Tc versus the spiral wavevector Q for different values 
of the interface parameter 7;,. Here df = 5 nm and d s = 7 
nm. The remaining parameters are the same as in Fig. 2 and 
3 of Ref. El 





FIG. 5: Singlet and triplet pair amplitudes at the S/F in- 
terface versus the spiral wavevector Q. An increase of Q is 
accompanied by an enhancement of the singlet amplitude to- 
gether with a reduction of the (negative) triplet amplitude. 
Here, we took the same parameters as in Fig.|3]with df = 2£ s . 



tcrfacc characterized here by 7& (in the case of a perfect 
transparency at the interface jj, = 0, while 3> 1 for 
a low-barrier transparency). From the Fig. 0] one sees 
that a relatively high transparency of the S/F interface 
is needed in order to have a strong dependence of T c on 
the wavevector Q. 

So far we have only discussed the evolution of T c with 
the exchange field inhomogencity. In the following we 
show how the singlet and triplet pair amplitudes depend 
on the spiral wavevector Q. For this purpose, we study 
the quantities 

Fs(z) = T22f s (u} n ,z) 

n>0 

Ftz(z) = Ty^^Imftzjun, z). 

n>0 

In the Fig. [3 one sees that the enhancement of T c with Q 
can be related with the enhancement of the singlet am- 
plitude J- s (0) at the interface. We considered the param- 
eter 7f, = 0, which means that the singlet and triplet pair 
amplitudes are continuous at the interface. The triplet 
amplitude Ttz{Q) penetrating the singlet superconductor 
(which is negative in the plot) decreases to zero when Q 
increases. Thus, it seems that a reduction of the triplet 
pair correlations is in correspondence with an enhance- 
ment of the superconducting critical temperature T c . 

The most important effect of the inhomogencity con- 
cerns the features of the spatial dependence of the pair 
amplitudes. It is known that a fixed exchange field may 
give rise to a change of sign of the pair amplitude in the 
ferromagnet at some distances of the S/F interface (this 
occurs under specific conditions, e.g. a sufficiently large 
df). Actually, both singlet and triplet amplitudes may 
have a non-monotonic spatial dependence and change 



FIG. 6: Spatial dependence of singlet and triplet amplitudes 
for Q^ 3 = and Q£ s = 1. Here we used the same parameters 
as in Fig. |3] with df = 2£ s . 



sign. This feature is illustrated in Fig. the curves 
for the singlet and triplet correlations change both sign 
in the ferromagnet at Q = 0. One can even note that 
they intersect near the position z = — 1.3 £ s . As shown 
in Fig. El the shape of the spatial dependences for the 
singlet and triplet amplitudes is modified in the presence 
of the inhomogencity (Q ^ 0). The crossing of the singlet 
and triplet curves tends to be suppressed with Q. Above 
a threshold value for Q, one finds that the singlet and 
triplet amplitudes do not change sign. 

Finally, it is possible to capture the numerical results 
concerning the reentrance of T c or the suppression of the 
nonmonotonic spatial dependence of the pair amplitudes 
with the exchange field inhomogeneity at a qualitative 
level. In fact, these features stem from the imaginary 
character of the eigenvalues k 2 which are for Q = 



2(u n ±iJ) /D f . 



(80) 



In Section III.B. we obtained that in the presence of the 
inhomogeneity, J has to be replaced by an effective ex- 
change field J±(Q) 



J±(Q) = JUl-rfTiri 



(81) 



where we remind that 77 = DfQ 2 /4J. It is worth noting 
that now the imaginary part of k 2 for 77 < 1 is propor- 
tional to ^/l — V 2 - Obviously, it is reduced in the pres- 
ence of the inhomogeneity of the exchange field (?/ ^ 0) 
compared to the value obtained for a fixed exchange field 
(i] = 0). Furthermore, for 77 > 1, the eigenvalues k^ ± 
become only real. This indicates that the inhomogeneity 
of the exchange field is detrimental to all the interesting 
features that characterize the S/F proximity effect such 
as the non-monotonic dependence T c (df), the oscillatory 
behavior of the pair amplitudes or of the local density of 
states in the ferromagnet. 
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IV. CONCLUSION 

We have discussed the system of spin-dependent qua- 
siclassical equations that describes diffusive S/F hybrid 
structures. We have pointed out that it is possible to 
predict many physical features of the S/F proximity ef- 
fect directly from these equations, independent of the 
specific geometry for the hybrid structure or the specific 
model for the spatial evolution of the exchange field in 
the ferromagnet. Singlet and triplet pair correlations are 
demonstrated to always coexist near the S/F interface. 
The triplet vector has the tendency to be aligned with 
the exchange field, resulting predominantly into triplet 
pairs with zero-spin projection. We have shown that the 
triplet correlations with nonzero spin projection on the 
local exchange field, which may have a long-range pene- 
tration in the ferromagnet, arise when perfect alignment 
between the triplet vector with the magnetic moment 
is prevented. The production of triplet components is 
accompanied by a spin-splitting in the local density of 
states and by the induction of a spin magnetization near 
the S/F interface. 

We have studied quantitatively the dependences of the 
singlet and triplet pair amplitudes and of the supercon- 



ducting critical temperature on the spiral wavevector 
within a model of a rotating magnetization in the fer- 
romagnet. The typical non monotonic behaviors of these 
quantities present for bilayers with homogeneously mag- 
netized ferromagnets are suppressed with increasing de- 
gree of inhomogeneity. This demonstrates the influence 
of domain walls in the ferromagnet on measurable prop- 
erties, implying that it is necessary to characterize the 
domain structure in order to compare quantitatively ex- 
periment with theory. The dependence of the supercon- 
ducting properties on the inhomogeneity of the exchange 
field is expected to be a general feature of the proximity 
effect in mesoscopic hybrid structures composed of super- 
conductors and ferromagnets, because it is a signature of 
the triplet correlations with zero-spin projection induced 
near the S/F interface. 
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